Telomere sequence variability in genotypes from natural plant populations: unusual block-organized double-monomer terminal telomeric arrays.

BACKGROUND: Telomeres are the nucleoprotein complexes that physically cap the ends of eukaryotic chromosomes. Most plants possess Arabidopsis-type telomere sequences (TSs). In addition to terminal TSs, more diverse interstitial TSs exists in plants. Although telomeres have been sufficiently studied, the actual diversity of TSs in land plants is underestimated. RESULTS: We investigate genotypes from seven natural populations with contrasting environments of four Chenopodium species to reveal the variability in TSs by analyzing Oxford Nanopore reads. Fluorescent in situ hybridization was used to localize telomeric repeats on chromosomes. We identified a number of derivative monomers that arise in part of both terminal and interstitial telomeric arrays of a single genotype. The former presents a case of block-organized double-monomer telomers, where blocks of Arabidopsis-type TTTAGGG motifs were interspersed with blocks of derivative TTTAAAA motifs. The latter is an integral part of the satellitome with transformations specific to the inactive genome fraction. CONCLUSIONS: We suggested two alternative models for the possible formation of derivative monomers from telomeric heptamer motifs of Arabidopsis-type. It was assumed that derivatization of TSs is a ubiquitous process in the plant genome but occurrence and frequencies of derivatives may be genotype-specific. We also propose that the formation of non-canonical arrays of TSs, especially at chromosomal termini, may be a source for genomic variability in nature.

A Multigene Phylogeny of Native American Hawkweeds (Hieracium Subgen. Chionoracium, Cichorieae, Asteraceae): Origin, Speciation Patterns, and Migration Routes.

Native American hawkweeds are mainly mountainous species that are distributed all over the New World. They are severely understudied with respect to their origin, colonization of the vast distribution area, and species relationships. Here, we attempt to reconstruct the evolutionary history of the group by applying seven molecular markers (plastid, nuclear ribosomal and low-copy genes). Phylogenetic analyses revealed that Chionoracium is a subgenus of the mainly Eurasian genus Hieracium, which originated from eastern European hawkweeds about 1.58-2.24 million years ago. Plastid DNA suggested a single origin of all Chionoracium species. They colonized the New World via Beringia and formed several distinct lineages in North America. Via one Central American lineage, the group colonized South America and radiated into more than a hundred species within about 0.8 million years, long after the closure of the Isthmus of Panama and the most recent uplift of the Andes. Despite some incongruences shown by different markers, most of them revealed the same crown groups of closely related taxa, which were, however, largely in conflict with traditional sectional classifications. We provide a basic framework for further elucidation of speciation patterns. A thorough taxonomic revision of Hieracium subgen. Chionoracium is recommended.

The structural diversity of CACTA transposons in genomes of Chenopodium (Amaranthaceae, Caryophyllales) species: specific traits and comparison with the similar elements of angiosperms.

BACKGROUND: CACTA transposable elements (TEs) comprise one of the most abundant superfamilies of Class 2 (cut-and-paste) transposons. Over recent decades, CACTA elements were widely identified in species from the plant, fungi, and animal kingdoms, but sufficiently studied in the genomes of only a few model species although non-model genomes can bring additional and valuable information. It primarily concerned the genomes of species belonging to clades in the base of large taxonomic groups whose genomes, to a certain extent, can preserve relict and/or possesses specific traits. Thus, we sought to investigate the genomes of Chenopodium (Amaranthaceae, Caryophyllales) species to unravel the structural variability of CACTA elements. Caryophyllales is a separate branch of Angiosperms and until recently the diversity of CACTA elements in this clade was unknown. RESULTS: Application of the short-read genome assembly algorithm followed by analysis of detected complete CACTA elements allowed for the determination of their structural diversity in the genomes of 22 Chenopodium album aggregate species. This approach yielded knowledge regarding: (i) the coexistence of two CACTA transposons subtypes in single genome; (ii) gaining of additional protein conserved domains within the coding sequence; (iii) the presence of captured gene fragments, including key genes for flower development; and (iv)) identification of captured satDNA arrays. Wide comparative database analysis revealed that identified events are scattered through Angiosperms in different proportions. CONCLUSIONS: Our study demonstrated that while preserving the basic element structure a wide range of coding and non-coding additions to CACTA transposons occur in the genomes of C. album aggregate species. Ability to relocate additions inside genome in combination with the proposed novel functional features of structural-different CACTA elements can impact evolutionary trajectory of the host genome.

A global survey reveals a divergent extradiol dioxygenase clade as a widespread complementary contributor to the biodegradation of mono- and polycyclic aromatic hydrocarbons.

Extradiol dioxygenation is a key reaction in the microbial aerobic degradation of mono- and polycyclic aromatic hydrocarbon catecholic derivatives. It has been reported that many bacterial enzymes exhibiting such converging functions act on a wide range of catecholic substrates. The present study reports a new subfamily of extradiol dioxygenases (EXDOs) with broad substrate specificity, the HrbC EXDOs. The new clade belongs to the XII cluster within family 2 of the vicinal oxygen chelate superfamily (EXDO-VC2), which is typically characterized by a preference for bicyclic substrates. Coding hrbC orthologs were isolated by activity-based screening of fosmid metagenomic libraries from large DNA fragments derived from heavily PAH-contaminated soils. They occurred as solitary genes within conserved sequences encoding enzymes for amino acid metabolism and were stably maintained in the chromosomes of the Betaproteobacteria lineages harboring them. Analysis of contaminated aquifers revealed coexpression of hrbC as a polycistronic mRNA component. The predicted open reading frames were verified by cloning and heterologous expression, confirming the expected molecular mass and meta-cleavage activity of the recombinant enzymes. Evolutionary analysis of the HrbC protein sequences grouped them into a discrete cluster of 1,2-dihydroxynaphthalene dioxygenases represented by a cultured PAH degrader, Rugosibacter aromaticivorans strain Ca6. The ecological importance and relevance of the new EXDO genes were confirmed by PCR-based mapping in different biogeographical localities contaminated with a variety of mono- and polycyclic aromatic compounds. The cosmopolitan distribution of hrbC in PAH-contaminated aquifers supports our hypothesis about its auxiliary role in the degradation of toxic catecholic intermediates, contributing to the composite EXDO catabolic capacity of the world's microbiomes.

Molecular Evolution and Organization of Ribosomal DNA in the Hawkweed Tribe Hieraciinae (Cichorieae, Asteraceae).

Molecular evolution of ribosomal DNA can be highly dynamic. Hundreds to thousands of copies in the genome are subject to concerted evolution, which homogenizes sequence variants to different degrees. If well homogenized, sequences are suitable for phylogeny reconstruction; if not, sequence polymorphism has to be handled appropriately. Here we investigate non-coding rDNA sequences (ITS/ETS, 5S-NTS) along with the chromosomal organization of their respective loci (45S and 5S rDNA) in diploids of the Hieraciinae. The subtribe consists of genera Hieracium, Pilosella, Andryala, and Hispidella and has a complex evolutionary history characterized by ancient intergeneric hybridization, allele sharing among species, and incomplete lineage sorting. Direct or cloned Sanger sequences and phased alleles derived from Illumina genome sequencing were subjected to phylogenetic analyses. Patterns of homogenization and tree topologies based on the three regions were compared. In contrast to most other plant groups, 5S-NTS sequences were generally better homogenized than ITS and ETS sequences. A novel case of ancient intergeneric hybridization between Hispidella and Hieracium was inferred, and some further incongruences between the trees were found, suggesting independent evolution of these regions. In some species, homogenization of ITS/ETS and 5S-NTS sequences proceeded in different directions although the 5S rDNA locus always occurred on the same chromosome with one 45S rDNA locus. The ancestral rDNA organization in the Hieraciinae comprised 4 loci of 45S rDNA in terminal positions and 2 loci of 5S rDNA in interstitial positions per diploid genome. In Hieracium, some deviations from this general pattern were found (3, 6, or 7 loci of 45S rDNA; three loci of 5S rDNA). Some of these deviations concerned intraspecific variation, and most of them occurred at the tips of the tree or independently in different lineages. This indicates that the organization of rDNA loci is more dynamic than the evolution of sequences contained in them and that locus number is therefore largely unsuitable to inform about species relationships in Hieracium. No consistent differences in the degree of sequence homogenization and the number of 45S rDNA loci were found, suggesting interlocus concerted evolution.

Characterization and Dynamics of Repeatomes in Closely Related Species of Hieracium (Asteraceae) and Their Synthetic and Apomictic Hybrids.

The repetitive content of the plant genome (repeatome) often represents its largest fraction and is frequently correlated with its size. Transposable elements (TEs), the main component of the repeatome, are an important driver in the genome diversification due to their fast-evolving nature. Hybridization and polyploidization events are hypothesized to induce massive bursts of TEs resulting, among other effects, in an increase of copy number and genome size. Little is known about the repeatome dynamics following hybridization and polyploidization in plants that reproduce by apomixis (asexual reproduction via seeds). To address this, we analyzed the repeatomes of two diploid parental species, Hieracium intybaceum and H. prenanthoides (sexual), their diploid F1 synthetic and their natural triploid hybrids (H. pallidiflorum and H. picroides, apomictic). Using low-coverage next-generation sequencing (NGS) and a graph-based clustering approach, we detected high overall similarity across all major repeatome categories between the parental species, despite their large phylogenetic distance. Medium and highly abundant repetitive elements comprise ∼70% of Hieracium genomes; most prevalent were Ty3/Gypsy chromovirus Tekay and Ty1/Copia Maximus-SIRE elements. No TE bursts were detected, neither in synthetic nor in natural hybrids, as TE abundance generally followed theoretical expectations based on parental genome dosage. Slight over- and under-representation of TE cluster abundances reflected individual differences in genome size. However, in comparative analyses, apomicts displayed an overabundance of pararetrovirus clusters not observed in synthetic hybrids. Substantial deviations were detected in rDNAs and satellite repeats, but these patterns were sample specific. rDNA and satellite repeats (three of them were newly developed as cytogenetic markers) were localized on chromosomes by fluorescence in situ hybridization (FISH). In a few cases, low-abundant repeats (5S rDNA and certain satellites) showed some discrepancy between NGS data and FISH results, which is due partly to the bias of low-coverage sequencing and partly to low amounts of the satellite repeats or their sequence divergence. Overall, satellite DNA (including rDNA) was markedly affected by hybridization, but independent of the ploidy or reproductive mode of the progeny, whereas bursts of TEs did not play an important role in the evolutionary history of Hieracium.

The major satellite DNA families of the diploid Chenopodium album aggregate species: Arguments for and against the "library hypothesis".

Satellite DNA (satDNA) is one of the major fractions of the eukaryotic nuclear genome. Highly variable satDNA is involved in various genome functions, and a clear link between satellites and phenotypes exists in a wide range of organisms. However, little is known about the origin and temporal dynamics of satDNA. The "library hypothesis" indicates that the rapid evolutionary changes experienced by satDNAs are mostly quantitative. Although this hypothesis has received some confirmation, a number of its aspects are still controversial. A recently developed next-generation sequencing (NGS) method allows the determination of the satDNA landscape and could shed light on unresolved issues. Here, we explore low-coverage NGS data to infer satDNA evolution in the phylogenetic context of the diploid species of the Chenopodium album aggregate. The application of the Illumina read assembly algorithm in combination with Oxford Nanopore sequencing and fluorescent in situ hybridization allowed the estimation of eight satDNA families within the studied group, six of which were newly described. The obtained set of satDNA families of different origins can be divided into several categories, namely group-specific, lineage-specific and species-specific. In the process of evolution, satDNA families can be transmitted vertically and can be eliminated over time. Moreover, transposable element-derived satDNA families may appear repeatedly in the satellitome, creating an illusion of family conservation. Thus, the obtained data refute the "library hypothesis", rather than confirming it, and in our opinion, it is more appropriate to speak about "the library of the mechanisms of origin".

Transposons and satellite DNA: on the origin of the major satellite DNA family in the Chenopodium genome.

Extensive and complex links exist between transposable elements (TEs) and satellite DNA (satDNA), which are the two largest fractions of eukaryotic genome. These relationships have a crucial effect on genome structure, function and evolution. Here, we report a novel case of mutual relationships between TEs and satDNA. In the genomes of Chenopodium s. str. species, the deletion derivatives of tnp2 conserved domain of the newly discovered CACTA-like TE Jozin are involved in generating monomers of the most abundant satDNA family of the Chenopodium satellitome. The analysis of the relative positions of satDNA and different TEs utilizing assembled Illumina reads revealed several associations between satDNA arrays and the transposases of putative CACTA-like elements when an ~ 40 bp fragment of tnp2 served as the start monomer of the satDNA array. The high degree of identity of the consensus satDNA monomers of the investigated species and the tnp2 fragment (from 82.1 to 94.9%) provides evidence of the genesis of CficCl-61-40 satDNA family monomers from analogous regions of their respective parental elements. The results were confirmed via molecular genetic methods and Oxford Nanopore sequencing. The discovered phenomenon leads to the continuous replenishment of species genomes with new identical satDNA monomers, which in turn may increase species satellitomes similarity.

Evolutionary history and genetic diversity of apomictic allopolyploids in Hieracium s.str.: morphological versus genomic features.

PREMISE: The origin of allopolyploids is believed to shape their evolutionary potential, ecology, and geographical ranges. Morphologically distinct apomictic types sharing the same parental species belong to the most challenging groups of polyploids. We evaluated the origins and variation of two triploid taxa (Hieracium pallidiflorum, H. picroides) presumably derived from the same diploid parental pair (H. intybaceum, H. prenanthoides). METHODS: We used a suite of approaches ranging from morphological, phylogenetic (three unlinked molecular markers), and cytogenetic analyses (in situ hybridization) to genome size screening and genome skimming. RESULTS: Genotyping proved the expected parentage of all analyzed accessions of H. pallidiflorum and H. picroides and revealed that nearly all of them originated independently. Genome sizes and genome dosage largely corresponded to morphology, whereas the maternal origin of the allopolyploids had no discernable effect. Polyploid accessions of both parental species usually contained genetic material from other species. Given the phylogenetic distance of the parents, their chromosomes appeared only weakly differentiated in genomic in situ hybridization (GISH), as well as in overall comparisons of the repetitive fraction of their genomes. Furthermore, the repeatome of a phylogenetically more closely related species (H. umbellatum) differed significantly more. CONCLUSIONS: We proved (1) multiple origins of hybridogeneous apomicts from the same diploid parental taxa, and (2) allopolyploid origins of polyploid accessions of the parental species. We also showed that the evolutionary dynamics of very fast evolving markers such as satellite DNA or transposable elements does not necessarily follow patterns of speciation.

Natural History of a Satellite DNA Family: From the Ancestral Genome Component to Species-Specific Sequences, Concerted and Non-Concerted Evolution.

Satellite DNA (satDNA) is the most variable fraction of the eukaryotic genome. Related species share a common ancestral satDNA library and changing of any library component in a particular lineage results in interspecific differences. Although the general developmental trend is clear, our knowledge of the origin and dynamics of satDNAs is still fragmentary. Here, we explore whole genome shotgun Illumina reads using the RepeatExplorer (RE) pipeline to infer satDNA family life stories in the genomes of Chenopodium species. The seven diploids studied represent separate lineages and provide an example of a species complex typical for angiosperms. Application of the RE pipeline allowed by similarity searches a determination of the satDNA family with a basic monomer of ~40 bp and to trace its transformation from the reconstructed ancestral to the species-specific sequences. As a result, three types of satDNA family evolutionary development were distinguished: (i) concerted evolution with mutation and recombination events; (ii) concerted evolution with a trend toward increased complexity and length of the satellite monomer; and (iii) non-concerted evolution, with low levels of homogenization and multidirectional trends. The third type is an example of entire repeatome transformation, thus producing a novel set of satDNA families, and genomes showing non-concerted evolution are proposed as a significant source for genomic diversity.

Decontamination by Persteril 36 may affect the reliability of DNA-based detection of biological warfare agents-short communication.

Persteril 36 is a disinfectant with a broad spectrum of antimicrobial activity. Because of its bactericidal, virucidal, fungicidal, and sporicidal effectiveness, it is used as a disinfectant against biological warfare agents in the emergency and army services. In case of an attack with potentially harmful biological agents, a person's gear or afflicted skin is sprayed with a diluted solution of Persteril 36 as a precaution. Subsequently, the remains of the biological agents are analyzed. However, the question remains concerning whether DNA can be successfully analyzed from Persteril 36-treated dead bacterial cells. Spore-forming Bacillus subtilis and Gram-negative Pseudomonas aeruginosa and Xanthomonas campestris were splattered on a camouflage suit and treated with 2 or 0.2 % Persteril 36. After the disinfectant vaporized, the bacterial DNA was extracted and quantified by real-time PCR. A sufficient amount of DNA was recovered for downstream analysis only in the case of spore-forming B. subtilis treated with a 0.2 % solution of Persteril 36. The bacterial DNA was almost completely destroyed in Gram-negative bacteria or after treatment with the more concentrated solution in B. subtilis. This phenomenon can lead to false-negative results during the identification of harmful microorganisms.

Functional adaptation of microbial communities from jet fuel-contaminated soil under bioremediation treatment: simulation of pollutant rebound.

To investigate the link between the functionality and the diversity of microbial communities under strong selective pressure from pollutants, two types of mesocosms that simulate natural attenuation and phytoremediation were generated using soil from a site highly contaminated with jet fuel and under air-sparging treatment. An increase in the petroleum hydrocarbon concentration from 4900 to 18,500 mg kg(-1) dw soil simulated a pollutant rebound (postremediation pollutant reversal due to residual contamination). Analysis of soil bacterial communities by denaturing gradient gel electrophoresis of PCR-amplified 16S rRNA gene fragments showed stronger changes and selection for a phylogenetically diverse microbial population in the mesocosms with pollutant-tolerant willow trees. Enumeration of the main subfamilies of catabolic genes characteristic to the site detected a rapid increase in the degradation potential of both systems. A marked increase in the abundance of genes encoding extradiol dioxygenases with a high affinity towards various catecholic substrates was found in the planted mesocosms. The observed adaptive response to the simulated pollutant rebound, characterized by increased catabolic gene abundance, but with different changes in the microbial structure, can be explained by functional redundancy in biodegrading microbial communities.

Metagenomics reveals diversity and abundance of meta-cleavage pathways in microbial communities from soil highly contaminated with jet fuel under air-sparging bioremediation.

The extradiol dioxygenase diversity of a site highly contaminated with aliphatic and aromatic hydrocarbons under air-sparging treatment was assessed by functional screening of a fosmid library in Escherichia coli with catechol as substrate. The 235 positive clones from inserts of DNA extracted from contaminated soil were equivalent to one extradiol dioxygenase-encoding gene per 3.6 Mb of DNA screened, indicating a strong selection for genes encoding this function. Three subfamilies were identified as being predominant, with 72, 55 and 43 fosmid inserts carrying genes, related to those encoding TbuE of Ralstonia pickettii PK01 (EXDO-D), IpbC of Pseudomonas sp. JR1 (EXDO-K2) or DbtC of Burkholderia sp. DBT1 (EXDO-Dbt), respectively, whereas genes encoding enzymes related to XylE of Pseudomonas putida mt-2 were not observed. Genes encoding oxygenases related to isopropylbenzene dioxygenases were usually colocalized with genes encoding EXDO-K2 dioxygenases. Functional analysis of representative proteins indicated a subcluster of EXDO-D proteins to show exceptional high affinity towards different catecholic substrates. Based on V(max)/K(m) specificity constants, a task-sharing between different extradiol dioxygenases in the community of the contaminated site can be supposed, attaining a complementary and community-balanced catalytic power against diverse catecholic derivatives, as necessary for effective degradation of mixtures of aromatics.